Spin-correlation functions in ultracold paired atomic-fermion systems: 
sum rules, self-consistent approximations, and mean fields 
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The spin response functions measured in multi-component fermion gases by means of rf transitions 
between hyperfine states are strongly constrained by the symmetry of the interatomic interactions. 
Such constraints are reflected in the spin f-sum rule that the response functions must obey. In 
particular, only if the effective interactions are not fully invariant in SU(2) spin space, are the 
response functions sensitive to mean field and pairing effects. We demonstrate, via a self-consistent 
calculation of the spin-spin correlation function within the framework of Hartree-Fock-BCS theory, 
how one can derive a correlation function explicitly obeying the f-sum rule. By contrast, simple one- 
loop approximations to the spin response functions do not satisfy the sum rule, except in special 
cases. As we show, the emergence of a second peak at higher frequency in the rf spectrum, as 
observed in a recent experiment in trapped 6 Li, can be understood as the contribution from the 
paired fermions, with a shift of the peak from the normal particle response proportional to the 
square of the BCS pairing gap. 

PACS numbers: PACS numbers: 03.75.Hh, 05.30.Jp, 67.40.Db, 67.40.Vs 
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I. INTRODUCTION 



Exploring pairing and superfluidity in ultracold 
trapped multicomponent-fermion systems poses consid- 
erable experimental and theoretical challenges [l], 0, H, 0, 
5]. Recently, Chin et al. 0] have found evidence, by rf 
excitation, of a pairing gap in a two-component trapped 
6 Li gas over a range of coupling strengths. The experi- 
ment, concentrating on the lowest lying hyperfine states, 
\<j) = |1), 1 2) and |3), with m F = 1/2, -1/2, and -3/2 
respectively, measures the long wavelength spin-spin cor- 
relation function, and is analogous to NMR experiments 
in superfluid 3 He @, H| ■ While at high temperatures the 
rf field absorption spectrum shows a single peak from 
unpaired atoms, at sufficiently low temperature a second 
higher frequency peak emerges, attributed to the contri- 
bution from BCS paired atoms. Theoretical calculations 
at the "one-loop" level of the spin response [l(| El 
support this interpretation. 

In this paper we carry out a fully self-consistent calcu- 
lation of the spin-spin correlation function relevant to the 
rf experiment, at the Hartree-Fock-BCS level, in order to 
understand the dependence of the response on mean field 
shifts and the pairing gap. The calculation requires going 
beyond the one-loop level, and summing bubbles to all 
orders, and is valid in the weakly interacting BCS regime, 
away from the BEC-BCS crossover - the unitarity limit. 
An important constraint on the mean field shifts was 
brought out by Leggett [l2j via a sum-rule argument: 
For a system with an interaction that is SU(2)-invariant 
in spin space, the spins in the long-wavelength limit sim- 
ply precess as a whole at the Larmor frequency, without 
mean field effects; then the spin-spin correlation function 
is dominated by a single pole at the Larmor frequency. 
While the effective interactions between the three lowest 
hyperfine states of 6 Li are not SU(2)-invariant the f-sum 
rule obeyed by the spin-spin correlation function still, as 



we shall show, implies strong constraints on the spin re- 
sponse, which are taken into account via a self-consistent 
calculation. 

In order to bring out the physics of a self-consistent 
approach to the spin response, we consider a spatially 
uniform system, and work within the framework of sim- 
ple BCS theory on the "BCS side" of the Feshbach reso- 
nance where the interactions between hyperfine states are 
attractive. We assume an effective Hamiltonian in terms 
of the three lowest hyperfine states explicitly involved in 
the experiments @, [Tj] (we take U = 1 throughout): 

H = f dr {T, (^ V ^( r ) ' + <^(r)Vv(r) 

1 3 

+ 2 E WV^(r)^(r)^(r)Vv(r)}, 

£J,fT / =l 

(1) 

where ipa- is the annihilation operator for state |cr), g aa > is 
the bare coupling constant between states a and a' , which 
we assume to be constant up to a cutoff A in momentum 
space. Consistent with the underlying symmetry we as- 
sume A to be the same for all channels, and take A — * oo 
at the end of calculating physical observables. The renor- 
malized coupling constants g aa i are related to those of the 
bare theory by 



9<7<7> ~ 9 aa ' 



mA/27r 2 



(2) 



where, in terms of the s-wave scattering length a aa > , 
g aa * = 47ra (T(T / /m. In evaluating frequency shifts in nor- 
mal states, we implicitly resum particle-particle ladders 
involving the bare couplings and generate the renormal- 
ized couplings. However, to treat pairing correlations 
requires working directly in terms of the bare g GG i [l4| . 

It is useful to regard the three states \a) as belonging 
to a pseudospin (denoted by Y) multiplet with the eigen- 
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values to ct of Y z equal to 1,0,-1 for a = 1,2,3. In terms of 
m a the Zeeman splitting of the three levels is 

4 = 4 - (4 - el)m„/2 + (4 + 4 - 2e%)m 2 j2. (3) 

The final term in 6 Li is of order 4% of the middle term 
on the BCS side. The interatomic interactions in the full 
Hamiltonian for the six F — 1/2 and F = 3/2 hyper- 
fine states are invariant under the SU(2) group of spin 
rotations generated by the total spin angular momentum 
F. The effective Hamiltonian can be derived from the 
full Hamiltonian by integrating out the upper three lev- 
els. However, because the effective interactions between 
the lower three levels depend on the non-SU(2) invariant 
coupling of the upper states to the magnetic field, the in- 
teractions in the effective Hamiltonian JT]) are no longer 
SU(2) invariant [HEl]. 

In the Chin et al. experiment equal numbers of 
atoms were loaded into states |1) and |2) leaving state 
|3) initially empty; transitions of atoms from |2) to 
|3) were subsequently induced by an rf field. Fi- 
nally the residue atoms in |2) were imaged, thus de- 
termining the number of atoms transferred to |3). 
The experiment (for an rf field applied along the 
x direction) basically measures the frequency depen- 
dence of the imaginary part of the correlation func- 

tion H)/d 3 r(r(^(r,i)V3(r, ^1(0,0)^(0,0))) (al- 

though in principle atoms can make transitions from |2) 
to 1 4); such a transition, at higher frequency, is beyond 
the range studied in the experiment, and is not of interest 
presently). Here T denotes time ordering. This correla- 
tion function can be written in terms of long-wavelength 
pseudospin-pseudospin correlation function, the Fourier 
transform of 



Xxx(t) = -i(T(Y x (t)Y x (0))); 



(4) 



here Y x = J d 3 ry x (r) is the x component of the total 
pseudospin of the system, 



2/x(r) = -^=(ipl(r)ip 2 (r) + ^4(r)t/>i(r) 
+V4(r)^ 3 (r)+^4(r)^ 2 (r)) 



(5) 



is the local pseudospin density along the x-axis. Since 
the experiment is done in a many-body state with N± — 
N2, the contribution from transitions between |1) and 
1 2) is zero 17]. The Fourier transform of Xxx(t) has the 
spectral representation, 



Xxx(ty 



7T 12 — U) 



(6) 



where Xxx( u ) = ^- m Xxx{u - «0+). 

In the next section we discuss the f-sum rule in gen- 
eral, review Leggett's argument, and illustrate how the 
sum rule works in simple cases. Then in Section III we 
carry out a systematic calculation, within Hartree-Fock- 
BCS theory, of the spin-spin correlation functions, gen- 
erating them from the single particle Green's functions. 



In addition to fulfilling the f-sum rule, our results are 
consistent with the emergence of the second absorption 
peak observed in the rf spectrum at low temperature from 
pairing of fermions. 



II. SUM RULES 

The f-sum rule obeyed by the pseudospin-pseudospin 
correlation function arises from the identity, 



+00 j, . 

—ux'M = m, H],Y X \). 



(7) 



The need for self-consistency is driven by the fact that 
the commutator on the right side eventually depends on 
the single particle Green's function, whereas the left side 
involves the correlation function. The static pseudospin 
susceptibility, — Xxx(0), is related to Xx X (uj) by 



Xxx(0) = 



d^x'Liu) 



(8) 



Leggett's argument that an SU(2) invariant system 
gives an rf signal only at the Larmor frequency is the 
following: Let us assume that the <? CTCT < are all equal, so 
that the interaction in Eq. ([1]) is SU(2) invariant in pseu- 
dospin space; in addition, let us assume, for the sake of 
the argument, that the Zeeman energy is —^fm a B z (7 
is the gyromagnetic ratio of the pseudospin). Then the 
right side of Eq. ((7]) becomes jB z (Y z ), while the static 
susceptibility, — Xxx{0), equals (Y z )/jB z . In this case, 
the spin equations of motion imply that the response is 
given by a single frequency (as essentially found experi- 
mentally @, HI). Thus for uj > 0, we take Xxxi^) to be 
proportional to 5(uj — ojq). Combining Eqs. ([7]) and ([5]), 
we find luq — "fB z , the Larmor frequency. The sum rule 
implies that neither mean field shifts nor pairing effects 
can enter the long wavelength rf spectrum of an SU(2) 
invariant system. 

It is instructive to see how the sum rule functions 
in relatively simple cases. We write the space and time 
dependent spin density correlation function as 

D xx (10) = -i(T (y x (l)y x (0))) 
= ±[£> 12 (1) + + £>2 3 (1) + An(l)], (9) 

where 

D Pa (l)= -i(T(^t(l)#(l)^( )^(0))), (10) 

and a, j3 = 1, 2, 3. Here -0(1), with 1 standing for {ri, t±}, 
is in the Hcisenberg representation, with Hamiltonian 
H' = H — J2 a ^<yN a . Equation |9]) implies that x"x i s a 
sum of x'pai wnere 

V 

X/3 a (^) = —D Pa (q = 7 n + fi a - fip), (11) 
and V is the system volume. 
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As a first example we consider free particles (denoted 
by superscript 0). For a ^ [3, 



D%{1) = -iG° a (-l)G°(l). 



(12) 



where G°(l), the free single particle Green's function, 
has Fourier transform, G° a (k, z) = l/(z — e£), with z the 
Matsubara frequency and e£ = k 2 /2m + e z — fi a . Then, 



1 N„-N, 



2 fi + e% - e 





z c z 



(13) 



from which we see that X /3ct( w ) ^ as a delta function peak 
at e§ — e z , as expected for free particles. This result is 
manifestly consistent with Eq. ([7]). 

Next we take interactions into account within the 
Hartree-Fock approximation (denoted by H) for the sin- 
gle particle Green's function, with an implicit resumma- 
tion of ladders to change bare into renormalizcd coupling 
constants. It is temp ting to factorize D as in the free 
particle case as d Efl EIES El, 

Df2(l) = -iG%(-l)Gf(l), 

where G% (k, z) = l/(z - Q), with 



(14) 



Cfe — o r c z 

2m 



(15) 



and n/3 the density of particles in hyperfine level (3. Then 

(16) 



n H0 (n - n o\ - n ° 7113 



and 

X/3a(0) = 
1 



(17) 



Consequently 



xjSaH = -^N fi -N a )S(u-AE Pa ), (18) 



where 



A^q, = e z + ^2 9f3a' n *' - e z - ^2 9a<jn a (19) 

is the energy difference of the single particle levels \a) 
and |/?). The response function x'pai^) is non-zero only 
at uj — AEp a . 

On the other hand, Xpai 00 ) obeys the sum rule 



= -j d 3 r([[^l(r)Mr),m^lm a (0)]) (20) 
= -{N ol -Np)[AEp a -g aP {ni3-n a )). (21) 



where the final line holds for the Hartree-Fock approxi- 
mation. The sum rule (|2ip is violated in this case unless 
g a /3 = 0. 

The self-consistent approximation for the correlation 
function (detailed in the following section) that maintains 
the sum rule and corresponds to the Hartree approxima- 
tion for the single particle Green's function includes a 
sum over bubbles in terms of the renormalized <?'s: 

dU*M = ,°™„£} or ( 22 ) 

1 +gp a D™(q,ll) 

Then with (l22l. 



X/J«» = \{N a - Np)5(u + (e z - e z 
+ 51 9a a n a - ^ gppUp +g a p(n/3 - n a )fj. (23) 

<r(#a) p(#/3) 

Note that x 3 2"M peaks at wj, = e 3 z -e% + (gi 3 - g 12 )n 1 , 
indicating that the mean field shift is (513 — 512)^1- 

This result agrees with the rf experiment done in a 
two level 6 Li system away from the resonance region [20| . 
This experiment finds that no matter whether the atoms 
in states |1) and \2) are coherent or incoherent, the rf 
signal of the transition between |1) and |2) never shows a 
mean field shift. As explained in [2Qj, in a coherent sam- 
ple, the internal degrees of freedom of all the fermions 
are the same, and thus there is no interaction between 
them. In the incoherent case, the above calculation gives 
Xi2"M = (^/2)(A r 2 - Ni)S(ui + e| - e z ), always peak- 
ing at the difference of the Zeeman energy, and therefore 
without a mean field contribution, plj 

In an rf experiment using all three lowest hyperfine 
states, the mean field shifts appear in X32" as (5i3 — 
gi2)nx. Since gaa' = 4:7ra aa >h 2 /m, our result (513—312)^1 
agrees with Eq. (1) of Ref. [la ]. However, from B =660 
to 900 G (essentially the region between the magnetic 
fields at which 013 and 012 diverge) no obvious deviation 
of the rf signal from the difference of the Zeeman energies 
is observed in the unpaired state [H, HI ■ The frequency 
shifts estimated from the result ([23]) taken literally in 
this region do not agree with experiment; one should 
not, however, trust the Hartree-Fock mean field approx- 
imation around the unitarity limit. The disappearance 
of the mean field shifts in the unitary regime has been 
attributed to the s-wave scattering process between any 
two different species of atoms becoming unitary-limited 
@, [I; however, the situation is complicated by the fact 
that the two two-particle channels do not become unitar- 
ity limited simultaneously. 



III. SELF-CONSISTENT APPROXIMATIONS 

References [11 and [11 laid out a general method to 
generate correlation functions self-consistently from the 
single particle Green's functions. To generate the corre- 
lation function Xxx{t), defined in Eq. J?]), we couple the 
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pseudospin to an auxiliary field F(r, t), analogous to the 
rf field used in the experiments, via the probe Hamilto- 
nian : 



iWbeW = -/ drF(r,t)y x {r). 



(24) 



The single particle Green's function is governed by the 
Hamiltonian H' , together with the probe Hamiltonian. 
The procedure is to start with an approximation to the 
single particle Green's function, and then generate the 
four-point correlation function by functional differentia- 
tion with respect to the probe field. Using this technique 
we explicitly calculate the pseudospin-pseudospin corre- 
lation functions consistent with the Hartree-Fock-BCS 
approximation for the single particle Green's function, in 
a three-component interacting fermion system, relevant 
to the rf experiment on the BCS side (g aC r> < 0). 

To calculate the right side of Eq. directly, we de- 
compose the Hamiltonian as H = H mvliI + H val , where 
-Hinvar is invariant under SU(2) and the remainder 



ff™ - 4 + (4 + 4 - 2e|)y//2 - ( e ; 

+ (<7l3 - 312) / V4VlV ; lV'3 + (923 - 912) 



4)^/2 

tp 3 W2 V>3, 
(25) 



is not invariant. We evaluate the right side of Eq. ||7J), 
([[Y x , H val ],Y x ]), term by term in the case that the states 
have particle number Ni — N2 = N and N 3 = 0, The 
Zeeman energy in H val gives (4 — 4)^> anc ^ the second 
term gives (512 - $13) /(V'lV'mV^)- 

We factorizc the correlation function within the 
Hartree-Fock-BCS theory for the contact pseudopoten- 
tial in ([1]), implicitly resumming ladders to renormalize 
the coupling constant in the direct and exchange terms 
[27l. to write, 



(512 - 513) J (^2 WlV>2) 



(312-313) /(^)(#) 

+(312-313) / (44)(i>ii>2). (26) 



Using Eq. @, we find a contribution from the second 
term, V(gi 2 - gi 3 )(n 2 ni + \A\ 2 / g 12 gi 3 ), where A = 
(4>i4>2)/gi2, is the BCS pairing gap between |1) and |2), 
assumed to be real and positive. The last term gives 

(312 - 323) / / /(V4 ^ 3 ^W 2 M 3 ^ 3 ^2) = °; altogether, 



7T 



(4 - e%)N - V{g 12 - g 13 ) (mn 2 + A 2 /g 12 g 13 ) . (27) 



gap contribute to the sum rule, allowing the possibility 
of detecting pairing via the rf absorption spectrum. 

We turn now to calculating the full Hartree-Fock-BCS 
pseudospin-pseudospin correlation function. For conve- 
nience we define the spinor operator 



* = (V'i,V'2,V'3,V'i,4>^ t 



3 ' 



and calculate the single particle Green's function 
G ab (l,l') = (-i){T* a (l)H>l(l')), 



(28) 



(29) 



where 1 denotes r\,t\, etc., and the subscripts a and b 
run from 1 to 6 (in the order from left to right in Eq. (128[) ; 
the subscripts 4, 5, and 6 should not be confused with the 
label for the upper three hyperfine states), and \& a (l) is 
in the Heisenberg representation with Hamiltonian H" = 
H' + H piohe (t). For F(r,t) = and with BCS pairing 
between |1) and |2), 



G = 



(G u 











G15 





\ 





G 22 





G24 














G33 
















G42 





G44 










G51 











G55 







V 














G66 


/ 



(30) 



To obtain a closed equation for G a t,(l,l'), we fac- 
torize the four-point correlation functions in the equa- 
tion of motion for G as before, treating the Hartree- 
Fock (normal propagator) and BCS (abnormal propa- 
gator) parts differently. In the dynamical equation for 
Gn (1,2), the term 0i 2 (V>| (1)^2 (1)^1(1)^4(2)) is approx- 
imated as gi2(V4(l)' ( / ; 2(l))('0i(l)' ! / ; i(2)) for the normal 
part, but g~i2 (ip2 (1)^1 (1))(V4 (1)V4(2)) for the abnormal 
part Q. Since ni = n 2 , e\ + g 12 n 2 + gi 3 n 3 - fi 2 = 
4 + 3i2"i + 323"-3 - M2 = — Mo> where fj, is the free 
particle Fermi energy; fj,Q enters into the single parti- 
cle Green's function as usual via the dispersion relation 

1/2 

Ek = [(fc 2 /2m — fi ) 2 + A 2 ] for the paired states. 

The equation of the single particle Green's function in 
matrix form is 



dl{G - 1 {ll)-F{l)TS(l-l) 

-£(11)}G(11') =5(1-1'), (31) 

where the inverse of the free single-particle Green's func- 
tion is 

GV^ll') = (iJL + |i ±A 5(1 - l')5 ab , (32) 



The absence of 1723 arises from N 3 being zero. Were all 
g acr ' equal, the right side of Eq. (|27]) would reduce to 
(4 — 4)^> as expected. When the interaction is not 

SU(2) invariant both mean field shifts and the pairing with the upper sign for a=l,2,3, and the lower for 
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=4,5,6. The matrix r is 



1 

71 



/O 1 








°\ 


1 1 











1 

















-1 








-1 





-1 


\o 





-1 


0/ 



The self energy takes the form 



(33) 



■iS(ll') 



/ 




312G12 


513G13 





— 312G15 


— 313G16 




P12G21 


— 912G11 — 523G33 


<723G23 


— 312G24 





— 323G26 




S13G31 


P23G32 


~ <?13Gn — (?23G22 


— <?13G34 


— 323G35 










512G51 


513G61 


512G22 + J?13G33 


— 312G21 


— 313G31 


V 


512 G42 





<723G62 


— 312G12 


gl2Gn + P23G33 


— 323G32 


1713 G43 


(723 G53 





— (713G13 


— 323G23 


513G11 + 323 



(34) 



where here G a b denotes G a f,(ll + )(5(1 — 1') with 1 + = 
We generate the correlation functions as 



D ab (12) = -n/2 



<5G afe (ll + ) 

^(2) J F=0 ' 



(35) 



(where the factor cancels that from the coupling of 
F(r, t) to the atoms via y x ) so that from Eq. (f3Tj) . 

+j^-(g,n))G(fc-g,z-n). (36) 



Using Eq. (f3"D"l) in (|3"6"|) . we derive 



^23 = 
D 12 = 



D 



2:i 



1 +.923^3 ^ 

^12 



! 



1 + 912^2 



where 



and 



D23 — II2233 + <7l3 



112433116651 
1 — 313II6611 



(37) 
(38) 

(39) 
(40) 



D 12 — II1122 — II1542; 
the bubble Tl a bcd(q, ^) is given by 

n Qbcd (q, fi)-^ G ab (k, z)G cd {k -q,z- fi), (41) 



and the summation on k is up to A. When A — > 0, 
Eqs. ((33) and ^ reduce to ([%?])■ since II2433 and n 6511 
are both proportional to A. Furthermore, when the in- 
teraction is SU(2) invariant, x'^i^) i s proportional to 
S(u> — (e^ — el)). If only 312 is non-zero, the response 
function D23 reduces to the single loop, II2233 (as cal- 
culated in Ref. and in fact satisfies the f-sum rule 
©■ 

To see that the result (13"T|) for the correlation function 
obeys the sum rule (f2"T)) . we expand Eq. (|3T|) as a power 
series in 1/57 in the limit 17 — > 00 and compare the coef- 
ficients of I/O 2 of both sides. In addition, with n\ — n 2 , 
we find J (cLl>/2'k)ujxi2{ uj ) = 0- 

Figure ([T]) shows the paired fermion contribution to 
X3 2 (o;), calculated from Eqs. ([37)1 and (jTTJ) , as a function 
of w, with g a .a' — ^kTi 1 a a ^i jm. This graph corresponds 
to the 6 Li experiment done in a spatially uniform sys- 
tem. The origin is the response frequency of unpaired 
atoms, which is uj^ 2 = e\ — e 2 . + (513 — 512)^1- We have 
not included the normal particle response in our calcula- 
tion and do not show this part of the response in the 
figure. The parameters used are kpa\ 2 = — tt/4 and 
fli3 = «23 = 0.1ai2, for which, T c = 0.084/iq. As the 
pairing gap, A, grows with decreasing temperature, the 
most probable frequency, w pa ir, in the response shifts to 
higher value. Within the framework of BCS theory, we 
can interpret the peak at higher frequency observed in 
the rf experiment as the contribution from the paired 
atoms. 

We now ask how the most probably frequency w pa i r is 
related to the pairing gap A. To do this we use the sum 
rule on x^fw), written in terms of X32( w )- Since 
Xxx(u) = X32M +X23H and X32 (-<*>) = -X2 3 M> we 
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-0.001 



0.003 



FIG. 1: (Color online). The pseudospin response function, 
1X32 1 s vs - u f° r fcFffli2 = — 7r/4 and (Z13 = 023 = 0.1ai2 for 
three temperatures. The curves correspond to a) T = 0.0831, 
A = 0.0084, b) T = 0.0830, A = 0.012, and c) T = 0.0820, 
A = 0.030, all in units of the free particle Fermi energy in 
hyperfme states |1) and |2). 



have 



f°° duj ,, . . n f 00 duj 11/ \ 

J-00 T J-oo 7T 



(42) 



Formally expanding Eq. ([57)) as a power series in I/O 
and comparing the coefficients of 1/17 on both sides, we 
find 



-xM = <n>/2. 

7T 



(43) 



Then assuming that the rf peak due to pairing is single 
and narrow (as found experimentally), we approximate 



X32M as tt(Y z )6(u> ■ 
pg]). we finally find 



r )/2. Using Eqs. fl27J), (gSJ) and 



w# = (gi3 —512)- 



A 2 



"0512513 



(44) 



where n = n% = n 2 . Thus BCS pairing shifts the spec- 
trum away from the normal particle peak by an amount 
proportional to A 2 . 

Equation (|44|) enables one to deduce the pairing gap A 
from the experimental data in the physical case, 313 7^ 0. 
However this result breaks down for the paired states 
when (713 = 0, a consequence of the dependence of 
the sum rule in Eq. (|27jl on the cutoff A of the bare 
model ([1]). To see this point, we note that the fac- 
tor (#13 — gxz) / giigxs that multiplies A 2 in Eq. pTj) 



arises from the combination of the bare coupling con- 
stants I/512 — 513/512! using Eq. jSJ we can write this 
combination in terms of the renormalized coupling con- 
stants and the cutoff as 



9l2 



toA 
2^ 



(g u L - mA/2^ 2 ) 



2\2 



(45) 

!l. :! »,.\ ■>-- 
Expanding in 1/A in the limit A — * 00, we see that the 
terms linear in A cancel, leaving the cutoff-independent 
result, — 5^, as in Eqs. (J27J and (|3I|). However, if 
only 512 is nonzero in this model then we hnd the cutoff- 
dependent result, 



A 2 A 2 , 1 

LU H = — = ( 

n gi2 n g 12 



raA, 
2^' 



(46) 



Fitting the measured shift in Ref. [6j, Fig. 2, to 
Eq. (|44|) . using the values of ai 2 and 013 as functions 
of magnetic field given in Ref. |13| . and assuming that 
g aa i = 4:nh 2 a crcr i /m, we find that for Fermi energy, 
E F = 3.6^K, A/E F = 0.23 at 904G {k F a 13 = -1.58, 
k F a±2 = —3.92), and 0.27 at 875G (or in terms of the 
Fermi momentum, k F ai3 = —1.69, k F cii2 — —6.31). 
Similarly for E F = 1.2^K, A/E F = 0.14 at 904G 
(k F a 13 = -0.91, fc F ai 2 = -2.26), and 0.19 at 875G 
(fc_Fai3 = —0.98, k F a\2 — —3.64). These values are in 
qualitative agreement with theoretical expectations [28| . 
although we expect corrections to the result (|44|) in the 
regime where the k F a are not small, and in finite trap 
geometry. 



IV. CONCLUSION 

As we have seen, the experimental rf result on the BCS 
side can be understood by means of a self-consistent cal- 
culation of the pseudospin response within the framework 
of BCS theory in the manifold of the lowest three hyper- 
fine states. The second peak observed at low temperature 
arises from pairing between fermions, with the displace- 
ment of the peak from the normal particle peak propor- 
tional to the square of the pairing gap A. The shift of the 
peak vanishes if the interactions within the lowest three 
states is SU(2) invariant. Although the results given here 
are for the particular case of the lowest hyperfme states 
in 6 Li, the present calculation can be readily extended 
to other multiple component fermion systems, as well as 
extended to include effects of the finite trap in realistic 
experiments. 
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